Remarks on the structure of Clifford quantum cellular automata 



Dirk-Michael Schlingemann 

ISI Foundation, Quantum Information Theory Group 
Viale S. Severn 65 
10133 Torino, Italy 

December 4, 2008 

Abstract. We report here on the structure of reversible quantum cellular automata with the additional 
restriction that these are also Clifford operations. This means that tensor products of Weyl operators 
(projective representation of a finite abelian symplectic group) are mapped to multiples of tensor 
products of Weyl operators. Therefore Clifford quantum cellular automata are induced by symplectic 
cellular automata in phase space. We characterize these symplectic cellular automata and find that all 
possible local rules must be, up to some global shift, reflection invariant with respect to the origin. In 
the one dimensional case we also find that all ID Clifford quantum cellular automata are generated by 
a few elementary operations. 



1. Introduction 

A standard modeling technique for various complex systems are cellular au- 
tomata. In fact they are ideally suited for models of diverse phenomena as coffee 
percolation, highway traffic and oil extraction from porous media. Cellular au- 
tomata also provide an abstract computational model that can simulate Turing ma- 
chines, and even explicit simple automata such as Conway's life game have been 
shown to support universal computation [2]. Quantum cellular automata provide 
model for analyzing quantum computational processes and quantum computational 
complexity. In his famous paper [3], Feynman discusses this idea in order to ob- 
tain a model for quantum computing which can be more powerful than a classical 
computer. Quantum cellular automata also play a role in view of quantum compu- 
tational complexity. This aspect has been studied by Bernstein and Vazirani [ 1 ] by 
using the concept of quantum Turing machines. Watrous (see e.g. [11]) continued 
this discussion by relating quantum Turing machines to quantum cellular automata. 
Last but not least, there may also be interesting applications besides quantum infor- 
mation theory, for instance, quantum cellular automata could serve as ultra-violet 
regularized quantum field theories. 

In this paper, we are concerned with the "quantized version" of cellular automata 
based on the concepts that are outlined in the article of Schumacher and Werner O. 



2 



D.-M. Schlingemann, Remarks on the structure of Clifford quantum cellular automata 



2. The general concept 

In order to motivate the concept for "quantum" cellular automata, we briefly 
review here the idea of reversible classical cellular automata from an algebraic point 
of view. In this context a complex classical system consists of single cells that are 
labeled by a countable set X. For many applications this set is given by a regular 
cubical lattice X = Z rf of dimension d. In our later discussion, we restrict our 
considerations to this case. To each cell x € X a finite set of states Q is associated, 
so that the classical configuration space of the total system is given by all functions 
from X into the set Q. A further important aspect is concerned with the local action 
of cellular automata which means that when the automaton is applied, the updated 
state of a single cell x e X only depends on the states of a finite set N(x) C X of 
"neighboring" cells. In the case of a regular cubic lattice X = Z d , the neighboring 
scheme is usually chosen to be translationally invariant. Here one takes some finite 
set N C 1j d and defines the neighbors of a cell x according to N(x) = N + x. 
In mathematical terms, a reversible classical quantum cellular automata is given 
a bijective map T: Q x — > Q x such that for each x € X there exists a function 
T x : Q N ^ -► Q which satisfies T(q)(x) = T x (q\ N{x) ). Here q is a function in Q 
and q\N(x) denotes the restriction of q to the neighborhood N(x). 

In the case of a regular lattice structure X = translation invariance is an 
additional requirement for the automaton T. The translation group Z d acts naturally 
on the configuration space according to (T x q)(y) = q(y — x). Translation 
invariance for T means that T commutes with all translations t x . As a consequence 
the cellular automaton T is completely determined by its "local rule" To: Q N — ► Q. 
Namely, the local rule T x at any cell x can be calculated from the local rule To at 
the origin by T x (q\ N+x ) = T (T- x q\ N ). 

In order to motivate the quantized concept, we reformulate the classical quan- 
tum cellular automata algebraically. If Q is a finite set, then the configuration space 
Q x is compact in the Tychonov topology. The observable algebra of the classical 
system is given by the abelian C*-algebra of continuous functions C(Q X ) which 
is canonically isomorphic to the tensor product C(Q X ) = ® x exC(Q), where 
C(Q) = is the abelian C*-algebra of functions on the single cell configuration 
space Q. In this picture, a reversible cellular automaton T induces an automorphism 
a on C(Q X ) by thepullback a(f)(q) = f(T(q)). 

The functions in C(Q X ) with values in the interval [0, 1] can be regarded as 
classical observables. Such an observable is localized in a subset U C X if the 
corresponding function / only depends on the restriction q\u of a classical config- 
uration q. With abuse of notation we express this fact as f(q) = f(q\u)- Thus if 
we restrict the automorphism a to observables that are localized in a single cell x 
the resulting observable is localized in the neighbor hood of x. This can be verified 
as follows: If an observable / is localized at x, then the value f(q) = f(q(x)) 
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only depends on q(x). The application of the automorphism a therefore gives 
«(/)(?(»)) = f{ T (<l){x)) = f(Tx(q\N(x))> where T x is the local rule at x. Thus 
the automorphism a propagates the localization region of an observable only into 
its neighboring cells. 

The basic idea to quantize the concept of cellular automata is to replace the clas- 
sical systems by quantum systems, i.e. theabelianC*-algebra(7((5 x ) = ®x^xC{Q) 
is replaced be a non-abelian one. The system under consideration is now described 
by a tensor product 



where to each cell x there is a (finite dimensional) C*-algebra 2l(x) associated with. 
We require here, that to each cell x an isomorphic copy of a fixed C*-algebra is 
assigned, i.e. %{x) = 2lo- Since X can be any countable set, we have to deal 
here with infinite tensor products. However, this is well defined in terms of the so 
called inductive limit. The algebra 21 is usually called the "quasi-local" algebra of 
observables. For each finite subset U C X there is a natural "local" subalgebra 
21(77) = ® x& u${(x) of the quasi local algebra. The operators a G 21(77) are identi- 
fied with operators in 2l(X) by filling the remaining tensor positions X\U with the 
unit operator. The concept of a reversible quantum cellular automaton is defined as 
follows: 

DEFINITION 1 . A reversible quantum cellular automata (QCA) is a *-automorphism 
a of the quasi local algebra 21 that fulfills the "locality condition": For each cell 
x G X and for each operator A G 2l(x) the operator a(A) G 2l(iV(x)) is localized 
in the neighborhood of x. 

As Definition [2] indicates, the concept of a QCA works for any type of lattice 
X with an appropriate neighborhood scheme where the locality requirement is the 
essential ingredient. In the subsequent analysis we focus on regular cubic lattices 
X = Z d only. In this case we have a natural action of the lattice translation group 
X d by automorphisms t x on the quasilocal algebra where t x is determined by 



with A x G 2l(x). We now consider those QCAs that respect the symmetry of lattice 
translations. 

DEFINITION 2. A translationally invariant reversible QCA is a reversible QCA a 
that commutes with the lattice translation group: a o t x = t x o a. 



21 = (g) 2l(x) 



(1) 
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In the subsequent we always refere to the translationally invariant situation. The 
translation symmetry can be exploited for the structural analysis of QCAs. In fact, 
a translationally invariant QCA is completely determined by its local rule at the 
origin. Recall that the local rule ao at x = is the restriction of the QCA a to the 
algebra 21(0). Due to the locality condition, there is a finite subset N C Z d such that 
a(2l(0)) C 2l(iV). To be compatible with translation invariance, the neighborhood 
scheme can be chosen such that N(x) = N + x and the "global rule", which is just 
the automorphism a, can be expressed in terms of the local rule ao by 




Thus a translationally invariant QCAs can be described in terms of its local rule 
only. In particular, if the single cell algebras are finite dimensional, the consruction 
of the QCA is a problem in finite dimensions. 

A strategy for constructing a QCA is based on finding a valid local rule. One has 
to choose a *4iomomorphism ao: 21(0) — > 2l(iV) and has to check the commutator 
condition [T x (a (A)) , a (B)] = for all A, B £ 21(0) and for all x G X with 
N n N + x ^ 0. Assuming that the single cell algebra 21(0) is finite dimensional, 
there are only finitely many conditions to be tested. For a comprehensive review on 
this issue, we refere here to the work of Schumacher and Werner J9). 

Although there are only finitely many conditions to check, a general systematic 
classification of QCAs is a highly non-trivial and still unsolved task. But there 
are particular classes of QCAs for which a complete and explicit classification is 
possible, as the class of Clifford (or quasifree) quantum cellular automata which we 
review here in the following. The results that we are presenting here are based on 
our previous article [8]. 

3. Clifford quantum cellular automata 

To explain the concept of Clifford quantum cellular, we consider a regular cubic 
lattice Jj d . To each cell we associate a full matrix algebra 2l(x) = M p (C) where p 
is a prime number. Moreover we choose a basis of Weyl operators in M p (C). These 
operators are generalizations of the Pauli operators and are constructed by shift and 
multiplier unitaries. To be more precise, we consider an orthonormal basis \ q) of the 
Hilbert space C p that is labled by elements q of the finite field ¥ p = 7L V . One way 
to define the Weyl operators is to determine its action on the basis | q) according to 

w(C)l?> = w(£+, £_)|?> = e p +q \q + £-> (4) 

where e p is the pth root of unity. As a consequence, the Weyl operators fulfill the 
relation 

w(£ + v) = 4-" + w(0wfa) (5) 
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which shows that the Weyl operators form a unitary projective representation of the 
additive group ¥ p . 

For the special case p = 2, which corresponds to qubits, the corresponding Weyl 
operators are related to the Pauli matrices X, Y, Z by X = w(0, 1), Z = w(l, 0) 
andy = -iw(l,l). 

We are now concerned with the quasi local algebra 21 which is given by the 
infinite tensor product of single cell algebras 2l(x) = M p (C) over the regular lattice 
Z d . The Weyl operators for the infinite system are given by tensor products 

w(0:=®w(£(a;)) (6) 

where the "phase space" vector £ is a function from the lattice Z d to the vector space 
Fp with finite support. Note that the finite support condition guarantees that only 
finitely many tensor factors are different from the identity which implies that w(£) 
is a well defined unitary operator that belongs to quasilocal algebra. Moreover, the 
complex linear hull of the Weyl operators is norm dense subalgebra of 21. I this 
sense, the Weyl operators form a "basis" for the quasi local algebra. 

Obviously, the Weyl operators of the infinite system fulfill the relation 

w(£ + r?) = e£&")w(0wfa) with /?(£,r?) = ^+(x)r ] -(x) . (7) 

which implies the commutation relation 

wfa)w(£) = w(£)wfa) with aft, n) = P(£, r,) - 0( v , . (8) 

Note that the symplectic form a for the infinite system is well defined since it is 
evaluated only for function with finite support. This relations justifies to interprete 
the functions £ as vectors in a discrete phase space — denoted by E p j in the fol- 
lowing — with symplectic form a. We are now prepared to give a precise definition 
of Clifford quantum cellular automata. 

DEFINITION 3. A Clifford quantum cellular automata (CQCA) a is a transla- 
tionally invariant reversible QCA which maps Weyl operators to multiples of Weyl 
operators. Thus there exists a function S: S Pj< j — > H p d as well as a phase- valued 
function ip: E p ^ h-> U(l) = {z G C\\z\ =1} such that 

a(w(0) = v(0w(S0 (9) 
holds for all phase space vectors £. 

A simple CQCA is given by lattice translations. The lattice translations act on 
the phase space vectors £ by just translating the function (r y £)(x) = £(x — y). With 
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abuse of notation we use the same symbol for the action on phase space as for the 
action on the quasilocal algebra. By construction the covariance relation 

t x (w(0) = w(r^) (10) 

follows immediately. Hence, the lattice translations t x are CQCAs in the sense of 
the definition given above. 

As we will see, the translation invariance together with the condition to map 
Weyl operators to multiple of Weyl operators is sufficient to determine a CQCA. 
This means, that the locality is a consequence of these conditions. To deal with the 
translation symmetry in an appropriate way, we identify the phase space Ep, d = 

d as a two dimensional module over the ring D p ^ of functions from the lattice 
Z d into the finite field ¥ p having finite support. The multiplication in the ring is the 
convolution of functions which is given by 

f*g = Y,f( x >*9- (ii) 

X 

We are now prepared to state the first structure theorem on CQCAs. We refere 
here the reader to our article 1 8 ] for a complete discussion of the proof in which 
uses techniques from the theory of projective representations of symplectic abelian 
groups (see e.g. lPT2ll ) as well as results from the theory of covariant completely 
positive maps |[T0ll5ll. 

THEOREM 4. For each CQCA a, there exists a two-by-two matrix s £ M2(D Ptq ) 
with entrees in the ring D p ^ and a translationally invariant phase valued function 
<p: D^ d —¥ U(l) such that 

a(w(£)) = p(£)w(s*£) (12) 

and if fulfills the cocycle condition 

<p(£ + v ) = e^-^ s ^<f(0<p(v) ■ (13) 

Moreover, the map s* preserves the symplectic form a, i.e. cr(s*£, s*7y) = rj(£, rj). 

We sketch here just the basic idea of the proof: It follows from the Weyl re- 
lations that each automorphism a that maps Weyl operator to multiples of Weyl 
operators according to (|9]> induces a F p -linear map S on phase space that preserves 
the symplectic form. Moreover, the condition to be an automorphism implies that 



the phase valued function ip fulfills ( 13 1 where, for this moment, we have to repace 
the operator s* by S. 

By taking advantage of the translation invariance, the phase valued function 
ip is translationally invariant and the symplectic map S commutes with the lattice 
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translations, it follows that S is given by the convolution s* with a matrix-valued 
function with finite support. Note that a two-by-two matrix s with entries in the 
ring D Pt d can equivalently be seen as a function that maps a lattice site x G Z d to a 
two-by-two matrix with entries in the finite field ¥ p . The convolution with a phase 
space £ vector is given by (s * £)(:c) = J2 y s (y) ' £( x ~ u) where • is usual matrix 
multiplication. The localization region of a Weyl operator w(£) is just the support 
of the function £. If the support of £ is just a single site x, then s * £ has support in 
N + x, where ./V is the support of s. Therefore, the application of the corresponding 
CQCA yields an operator a(w(£)) = <£>(£)w(s*£) which is localized in N + x. As 
a consequence, the support of s determines the neighborhood scheme of the QCA. 

We also shown in [8] that each matrix-valued function s for which the convo- 
lution s* preserves the symplectic form a a phase valued function 93 can be found, 
such that the equation ( p~2| ) defines a CQCA. Therefore, the classification of CQCAs 
is equivalent to characterize all two-by-two matrices s G M2(D P: d) with entries in 
the ring D p> d whose convolution s* is symplectic. In accordance with [§1, we call 
the convolution s* a "symplectic cellular automaton (SCA)". 

4. On the structure of Clifford quantum cellular automata 

For the further analysis of CQCAs, we have a closer look at the ring and module 
structure of the phase space S Pj ^ = d . As already mentioned E PjC i is a two 
dimensional module over the ring D p a where the product is the convolution. A 
function / G D Pi d acts on a phase space vector by / * £ = / * (£+,£-) = (/ * 
£ + , /*£_). A symplectic cellular automaton (SCA), which induces a CQCA, is then 
a module homomorphism. Namely, since the convolution in D P: d is commutative, 
we observe that s*/*£ = /*s*£. 

To take advantage of the translation symmetry in an appropriate manner, we 
have introduced the "algebraic Fourier transform", which identifies the ring D P) d 
with the commutative ring of Laurent polynomials 



generated by the variables u%, ■ ■ ■ ,u<i and its inverses u 1 , ■ ■ ■ , u . . For a function 
/ S D p ^ the corresponding Laurent-polynomial is simply given by 



where we write u x = u^u^ 2 ■ • • u x d d for a handy notation. To view the elements in 
D P; d as formalpolynomials, gives us a convenient book-keeping at hand. Namely, 
the convolution turns into a product of polynomials, i.e. for two functions the iden- 
tity 



D p ,d = ¥ p [ui,u 2 , ■ ■ ■ ,Ud,u x 



(14) 




(15) 



x 



f*h = fh 



(16) 
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holds. We mention here that the ring D p ^ is a "divison ring" which means that 
fh = implies that either / = or h = 0. Morover, the only invertible elements 
in D p ,d are the monomials u x with x G Z d . 

With help of this ring isomorphism /•—>/, the phase space can be identified 
with Dp d and, since a SCA is a module homomorphism, its Fourier transform just 
acts by matrix multiplication. To be more precise, a SQCA is given by a two-by-two 
matrix s with entries in the polynom ring D p ^ acting on a phase space vector by 

After applying the algebraic Fourier transform, the symplectic form is a -De- 
valued bilinear from E on Di , that can be calculated according to 

m,r])=Uv~-^r ]+ , (18) 

where / i— > / is the involution on the ring D p ^ which replaces in the polynomial 
/ = /(«) the variable by its inverse u^ 1 . In the lattice space, this corresponds to 
a reflection at the origin. The form E is related to the underlying symplectic form a 
by 

E(e,r ? ) = ^a(|,r^K (19) 

X 

where / i— > / is the inverse algebraic Fourier transform sending a polynomial / to 
a function / on the lattice. It is not difficult to observe, that E is a module homo- 
morphism in the second argument, i.e. E(£, frj) = E(£, rj)f and that it fulfills the 
relation E(£,f?) = — E(t/,£) = — E(^, £). Thus E is antisymmetric for reflection 
invariant polynomials. In this context, a helpful lemma for the characterization of 
SCAs is the following: 

LEMMA 5. A two-by-two matrix s G M2(-Dp^) with entries in the polynom ring 
D p d is a symplectic cellular automaton, if and only if, it preserves the form E, i.e. 
the identity E(s£, srj) = E(£, 77) holds. 

According to this lemma, the characterization of CQCAs (hence SCAs) reduces 
to the problem of finding two-by-two matrices with entries in the ring D p a which 
preserve E. Since we deal here with module homomorphism — a linear structure 
over the ring D p ^ — we have reduced a problem in infinitely many degrees of 
freedom to an effectively two-dimensional problem. 

There is an important subring in D p ^, which we denote here by P p ^,which 
consists of all polynomials that are invariant under the reflection Uk ► u^ 1 which 
means / = f(u) = f(u^ r ) = f. For d = 1, the corresponding function / in lattice 
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space is then given by a palindrome string / = (q n q n -i ■ ■ ■ qiqoqiq2 • • • q-n) starting 
at the left boundary of the support x = — n and ending at the right boundary x = n. 
For this reason, we call P PjC i the polynom subring of palindromes in D Pj d- If we look 
at the properties of the form, S, we see that it is a non-degenerate antisymmetric 
P p ^-bilinear form on P^ d . In analogy, that the symplectic group of C 2 is given by 
the special linear group SL(2, C) a first guess is, that the group of SCA is given 
by all two-by-two matrices with entries in the palindrome subring P p & having ring- 
valued determinant equal to one. Indeed, if we choose a matrix s € SL(2, Ppd), 
then we observe by a straight forward calculation that s preserves the form E. If 
we multiply s with a monomial u a , which corresponds to a lattice translation by 
a 6 Z rf , then we observe that £(it a s£, u a sr]) = £(s£, sr])u~ a u a = £(£, 77). Thus 
if s is a SCA then u a s is a SCA too. Indeed, all SCAs are of this type. The precise 
statement which we have established in ||8] Theorem 3.4] is the following: 

THEOREM 6. The group of Clifford quantum cellular automata acting on a d- 
dimensional lattice with single cell algebras M p (C) is isomorphic to the direct 
product Z d x SL(2, Ppd) of the lattice translation group and the special linear 
group of two-by-two matrices with entries in in the palindrome subring P P) d- 

This theorem can be used to build up a simple cooking recipe for constructing 
CQCAs. Firstly, take two palindromes /, h £ P P: d- Recall that palindromes are easy 
to get. Namely, for the case that q is not a palindrome you just make one by taking 
h = 9 + 9- Secondly, factorize the polynomial 1 — fg = f'h' in the subring P p ^ 
into two palindromes /', h! . Finally, you get your CQCA by building the matrix 



For this type of recipe, to find all possible factorizations of the polynomial 1 — fh 
is the crucial problem which can be quite cumbersome. However, there is allways 
the trivial solution which is given by h! = 1 and /' = 1 — fh. 

Concerning the factorization problem, at least for a one dimensional lattice d = 
1 the situation can be tackled. Here one takes advantage of the fact that the ring Pi p 
is a so called "Euclidean ring" (see e.g. 0) and an extended Euclidean algorithm 
for finding greatest common divisors can be applied. This yields in an even stronger 
classification result of one-dimensional CQCAs than provided by the dimension 
independent Theorem [4?] We have shown the following (8J Theorem 3.11]: 

THEOREM 7. Every Clifford quantum cellular automata s acting on a one dimen- 
sional lattice with single cell algebras M P (C) can be factorized into a product of a 
unique shift u a and elementary CQCAs of the following two types: The first type is 
a shear transformation 




(20) 




(21) 
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depending on an integer n € N. The second type depends on a constant c £ F p 
according to 



The automata f c are constant matrix valued polynomials which implies that they 
act on each lattice site independently. The have zero propagation speed since their 
neighborhood scheme only consits of the origin. On the other hand the propagation 
of the localization region of an observable is induced by the shear automata g n . 
Their local rule propagate from the origin into the cells {— n, n}. 



In this note, we have discussed some aspects on the structure of quantum cellular 
automata where we have mainly focused on our results on Clifford quantum cellular 
automata [El. 

We have characterized the group of CQCAs in terms of symplectic cellular 
automata on a suitable phase space. With the help of the concept of algebraic 
Fourier transform, this phase space can be identified with two-dimensional vec- 
tors of Laurent-polynomials, and symplectic cellular automata can be described by 
two-by-two matrices with Laurent-polynomial entries. We have reported that these 
entries must be reflection invariant and that up to some global shift the determi- 
nant of the matrix must be one, so the group of CQCAs is isomorphic to the direct 
product of the lattice translation group with the special linear group of two-by-two 
matrices with reflection invariant polynomials as matrix elements. 

Due to the specialty that for a ID lattice we are faced with an Euclidean ring, 
each one-dimensional CQCA can be factorized into a product of elementary shear 
automata and local transforms. 

Besides the core results, that we have presented here, there is a correspondence 
between ID CQCAs and ID translationally invariant stabilizer (graph) states (see 
e.g. flUEl for the notion of stabilizer (graph) states). For a fixed translationally in- 
variant pure stabilizer state, which is in particular a product state, every other trans- 
lationally invariant pure stabilizer state can be created by applying an appropriate 
CQCA. 

A further natural question is concerned with lattices with periodic boundary 
conditions. Here the techniques from infinitely extended lattices can be applied to a 
certain extend. The technical problem is here, however, that the involved polynom 
ring is no longer a division ring. 




(22) 



5. Concluding remarks 
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